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We calculate, as a function of temperature and density, the electron-electron interaction induced 
quasiparticle effective mass renormalization in 2D electron systems within the leading-order dy- 
namically screened Coulomb interaction expansion. We find an unexpected nonanalyticity and 
nonmonotonicity in the temperature dependent effective mass with the renormalized mass linearly 
increasing with temperature at low temperatures for all densities. 
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I. INTRODUCTION 

A key insight of the Landau Fermi Liquid theory is 
that interactions in a Fermi system lead to the renormal- 
ization of the single particle fermion mass giving rise to 
"quasiparticles" with renormalized effective mass whose 
low energy behavior is qualitatively similar to the corre- 
sponding noninteracting free particles^. Then, various 
single particles properties, e.g. specific heat, density of 
states, etc., in the interacting Fermi system are simply 
given, at least in the leading-order theory, by replacing 
the bare (i.e. "free particle") mass m by the correspond- 
ing renormalized effective mass m* . In this paper we 
present the very first microscopic calculation of the tem- 
perature dependent effective mass renormalization in an 
interacting 2D electron system (2DES), finding in the 
process an unexpected nonanalytic and nonmonotonic 
behavior of the effective mass m*(T) as a function of 
temperature in the 2DES. In particular, m*(T) first in- 
creases linearly with temperature in a 2DES reaching a 
density dependent maximum around T/Tp < 0.1 — 0.5, 
where Tp is the noninteracting Fermi temperature, af- 
ter which it decreases with increasing temperature. This 
nonmonotonic behavior, in particular the temperature 
induced enhancement of the 2DES quasiparticle effec- 
tive mass at low temperatures, is entirely unexpected be- 
cause the naive expectation is that quantum many-body 
electron-electron interaction effects (underlying the effec- 
tive mass renormalization phenomenon) should decrease 
with increasing temperature since the high temperature 
system is necessarily a classical system. The nonanalytic 
linear-T dependence of m*(T) is also quite unexpected 
since the usual fermionic Sommcrfeld thermal expansion 
always results in a quadratic temperature correction. 

Our work is partially motivated by the great deal of 
recent activity in semiconductor-based 2DES, e.g. Si in- 
version layers, GaAs heterostructures and quantum wells, 
etc. where the 2D carrier density can be varied (by tun- 
ing an external gate voltage), modifying the strength of 
the electron-electron interaction usually measured^ by 
the dimcnsionless parameter r s — me 2 / (Y?^Jim) with n 
being the 2D carrier density and m the bare (i.e. band) 
mass. The r s -parameteri in 3D metals (defined with 
respect to 3D densities) is typically 3 — 5 whereas in 



semiconductor 2DES r s could vary from 1 (or less) to 
20 (or higher), depending on the specific semiconductor 
system and carrier density being studied. Since the effec- 
tive mass renormalization scales with r s (small and large 
r s respectively corresponding to weakly and strongly in- 
teracting electron systems), one expects interesting and 
important many-body quasiparticle renormalization in 
2DES, particularly at large r s . It is therefore not sur- 
prising that the issue of the effective mass renormaliza- 
tion in 2DES has been extensively studied, both exper- 
imentally and theoreticallji^iSiL2iS over the last thirty 
years. All these theoretical studies of quasiparticle mass 
renormalization have, however, been restricted to T = 
both in the 2D 5 i 6 i 7 i 8 i 9 and 3D±°i system. While this zero- 
temperature restriction makes perfect sense in 3D sys- 
tems where the relevant Fermi temperature Tp = Ep /ks 
(defining the temperature scale for the electron system) 
is extremely high (Tp ~ 10 4 K in metals), it makes little 
sense for extremely low density 2DES of current inter- 
est n i 12 i 13 where Tp < IK, making T/T F ~ 1 in the 
experimental temperature range. The temperature (and 
density) dependent effective mass renormalization calcu- 
lation presented in this paper therefore takes on addi- 
tional significance because a number of recent experi- 
ments have reported large 2D effective mass renormal- 
izatio n 11 : 12 ! 13 at low densities and low temperatures. We 
note in this context that the 2D effective mass renor- 
malization m*/m in our finite temperature many-body 
theory is a function of two dimensionless parameter r s 
(oc n -1 / 2 ) and T/Tp (oc n^ 1 , since ksTp — nh 2 n/m in 
2DES), which are however not completely independent 
of each other (since they both depend on the electron 



dcnsity)-in particular, T/Tp 
ture and changing density. 



for a fixed tempera- 



The structure of this paper goes as following: In sec- 
tion^we present the theory for our effective mass calcu- 
lation. In section ITTT1 we provide our numerical results of 
our calculated effective mass as a function of r s and T. In 
section llVl we present the analytical results for effective 
mass in the r s Cl and T/Tp <C 1 limit. We conclude in 
section Ivl with a brief discussion. 
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II. THEORY 
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FIG. 1: The Feynman diagram for the self-energy. The cir- 
cles are polarization bubbles, the dashed lines the Coulomb 
interaction, and the solid lines the electron Green's function. 

We consider a 2DES interacting via the long range 
Coulomb interaction. The effective mass renormaliza- 
tion is microscopically calculated- from the electron self- 
energy function £(k, ivi) defined at the Matsubara imag- 
inary frequency iv\ and 2D wavevector k. To calculate 
the electron self-energy, we make the well-known "GW" 
approximation^* 1 ^* of a leading order expansion in the 
dynamically screened Coulomb interaction (the corre- 
sponding Feynman diagram for the self-energy is shown 
in Fig.^), obtaining (h = 1 throughout): 
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where V q = 2ire 2 /q is the 2D bare Coulomb potential, 
ii>x = i(2l + VjirksT and ioj n = i2nirkBT are the usual 
fermion/boson odd/even Matsubara frequencies (I, n in- 
tegers), £k = k 2 /(2m) — /i, fi the chemical potential, and 
e(k, iut n ) is the RPA dynamical dielectric function, given 
by the sum of the polarization bubble diagrams: 



e(k,iuj n ) = 1- V„Il(\s.,iu; n ), 



(2) 



with II(k, ioj n ) the electronic 2D polarizability. Within 
RPA, we have 



n(k,MJ„) = 2 



d 2 q n f (£ q ) - n F (£ q _ k ) 
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(3) 



where np(x) = l/(e x / T + 1) is the Fermi distribution 
function. 

The form of retarded polarizability II(k, lo) 
1 1 : k . iui n — > u!+iO + ) has been provided by previous work 4 
at zero temperature: 
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where fi is the chemical potential. The finite temperature 
form of retarded polarizability can obtained from Eq. Q 
by 



oo 

n(k,w, W T) = J dp! 



4Tcosh 2 (4^i) 



(5) 



The quasiparticle energy E k is obtained from the 
Dyson equation using the analytically continued retarded 
self-energy S(k, iv\ — > uj + i0 + ) = E(k, ui): 



E k = & + Re£(k, E k ) 



(6) 



Eq.@ is exact, while our GW-approximation is the 
first order perturbation expansion in the dynamically 
screened interaction. There has been much discussion 
on whether one should use exact Eq. JBJ for the effec- 
tive mass or the so-called on-shell approximation, keep- 
ing only first order interaction terms by taking the first 
order iteration of Eq. ©: 



£ k = 6c + Re£(k,&). 



(7) 



The on-shell approximation is expected to be more ac- 
curate within the GW scheme as it effectively accounts 
for some vertex corrections and obeys the Ward identi- 
ties. This approach has previously been used in 2D* and 
3Di* zero-temperature effective mass calculations, and is 
regarded to be better than solving the full Dyson equa- 
tion (Eq. 10 above). The two approaches are identical 
in the high-density limit r, < 1. For r s > 1, they give 
qualitatively similar but quantitatively different results. 
The effective mass can then be derived from the relation 
1/m* = (k~ x dEk/ dk)\ k _ k ^, remembering that the bare 

band mass m is given by 1/m = (k^d^k/dk) \ k _ k '■ 
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k — kr- 



where kp is the Fermi momentum for the non-interacting 
2DES. 

We use three different techniques in calculating the 
self-energy: frequency sum, frequency integration, and 
plasmon-pole approximation. The first two techniques 
are equivalent to each other, and correspond to different 
ways of doing the analytic continuation of the imaginary 
frequency self-energy. The frequency sum technique is 
explained in Refi^*, and the frequency integration tech- 
nique, also called spectral representation, in Ref»2. In the 
frequency sum method, the retarded self-energy is given 
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by 
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where jig(i) = 1 / (e^/ 7 " — 1 ) is the Bose distribution func- 
tion. For the frequency integration method, the retarded 
self-energy is 



Re£(k,w) 



d 2 q 

d 2 q _ 
(2n) 2 J 2n* Vq 

Ub{v) +n F (Cq-k 
V - (£ q _k - Oj) 



2tt g e(g, v) 



(10) 



The plasmon-pole approximation (PPA) is a simpler 
techniquaSiiSiii for carrying out the frequency sum in the 
RPA self-energy calculation by using a spectral pole (i.e. 
a delta function) ansatz for the dynamical dielectric func- 
tion e(k, lu): 

Imc- 1 (k,w) = Cfc[«(w-w fc )-«(w + w fc )]/2, (11) 

where the spectral weight Ck and the pole uik of the 
PPA propagator in Eq. Ijllfl are determined by using the 
Kramers-Kronig relation (i.e. causality) and the /-sum 
rule (i.e. current conservation). We mention that u>k in 
Eq. (|llf> does not correspond to the real plasmon dis- 
persion in the 2DES, but simulates the whole excitation 
spectra of the system behaving as an effective plasmon 
at low momentum and as the single-particle electron- 
hole excitation at large momentum, as constrained by 
the Kramers-Kronig relation and the /-sum rule. Details 
on the PPA are available in literatur e - 6 ! 16 , including the 
finite-temperature generalization 17 . The PPA, which is 
known 6 i 16 i 17 to give results close to the full RPA calcu- 
lation of self-energy, allows a trivial carrying out of the 
frequency sum in the retarded self-energy function lead- 
ing to: 
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We calculate the self-energy by carrying out the 2D 
momentum integration (Eq. |^1 E3 H2|> as well as the 



frequency sum (Eq. |5J) and the frequency integral (Eq. 
I10[l in order to obtain the quasiparticle effective mass 
(Eq. |HJ). We emphasize that our reason for carrying out 
our calculation of the electron self-energy by three dif- 
ferent techniques (RPA frequency sum and integration, 
and PPA) is to completely ensure the numerical accu- 
racy of the calculated temperature dependent effective 
mass by comparing the consistency among the three sets 
of results. This is particularly significant since there is 
no existing temperature-dependent effective mass calcu- 
lation in the literature for us to compare with. The fact 
that our three sets of results are consistent with each 
other (and we reproduce the existing^S* 7 ^^ T — ef- 
fective mass results from our finite temperature theory) 
provides compelling support for our conclusions in this 
paper. Since our results obtained in the three techniques 
are in good agreement, we will only show here our effec- 
tive mass results using RPA frequency sum method for 
the sake of brevity. 



III. NUMERICAL RESULTS 
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FIG. 2: Calculated T=0 effective mass as a function of r s in 
a 2DES in the high r 3 region. 

First, we present our extreme low temperature 2D re- 
sult (T/T F « 1(T 4 ) in Fig. Eland Fig.Gfl to be compared 
with the existing T = 2D resulta 5 ' 6 ! 7 ! 8 ' 9 , for m*(r s ) in 
the r s = — 10 range, showing that the effective mass 
renormalization could be almost as large as 5 for dilute 
r s ~ 10 2DES. We emphasize that the results presented 
in Fig. El and Fig. |3 based on the T — > limit of our 
finite temperature theory are in quantitative agreement 
with the existing T = 2D RPA effective mass calcu- 
lations 5 , which were, however, restricted to the r s (< 5) 
regime. 

In Fig. 21 we show our calculated 2D m*(T) as a func- 
tion of T/Tp for different values of the 2D interaction 
parameter r s (= 1 — 10). Fig. shows the effective mass 
temperature dependence at high densities. In the low 
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temperature region the effective mass first rises to some 
maximum, and then decreases as temperature increases. 
This peculiar behavior is present at all densities. The 
nonmonotonic trend is systematic, and the value of T/Tp 
where the effective mass peaks increases with increasing 
r s . The initial increase of m*(T) is linear in T/Tp as 

T — > 0, and the slope ^pi^r) is almost independent of r s 
for very small r s (< 1), but increases with r s for larger r s 
values. We mention that we get somewhat stronger tem- 
perature dependence (i.e. larger dm*/dT) in our PPA 
calculation (not shown here). 
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FIG. 4: Calculated 2D effective mass as a function of T/Tf 
for different high r s values: r s = 10 —* 1 from top to bottom. 
Note that Tf oc rj 2 , making the absolute temperature scale 
lower for higher r s values. 

In Fig.EJwe show the dependence of the effective mass 
renormalization as a function of the interaction parame- 
ter r s for a few values of fixed temperature (rather than 
fixed T/Tp, remembering that Tp oc rj 2 since Tp oc n 
and r s oc rT 1 ! 2 ). The calculated m*(r s ) for fixed T val- 



ues is quite striking: For low fixed values of T, m*/m 
initially increases with r s even faster than the corre- 
sponding T — result, eventually decreasing with r s 
at large enough values (where the corresponding T/Tp 
values become large enough). This nonmonotonic be- 
havior of m*(r s ) as a function of r s for fixed tempera- 
tures showing a temperature-dependent maximum (with 
the value of r s at which the m* peak occurs decreas- 
ing with increasing T as in Fig. EJ is complementary to 
the nonmonotonicity of m*(T) in Fig. 0] as a function of 
T/Tp (at fixed r s ) and arises from the relationship be- 
tween the dimensionless variables T/Tp (oc r~ 2 ) and r s 

— 1/2 

(oc T F ) due to their dependence on the carrier density 
(i.e. Tp oc n and r s oc n^ 1 / 2 ). 

One immediate consequence of our results shown in 
Figs. Hand El is that m*{T/T F ,r s ) = m*(T,n) in 2DES 
could show a strong enhancement at low (but finite) tem- 
peratures and low electron densities (large r s ). Compar- 
ing with the actual system parameters for 2D electrons in 
Si inversion layersii*i£ and GaAs heterostructuresi^ (and 
taking into account the quasi-2D form factor effects^ ne- 
glected in our strictly 2D calculation) we find that, con- 
sistent with recent experimental findingsiii^il^, our the- 
oretical calculations predict (according to Figs. 01 and El 
as modified by subband form factors) m*/m to be en- 
hanced by a factor of 2 — 4 for the experimental densities 
and temperatures used in recent measurementsiiii2ii£. 
Due to the approximate nature of our theory we do not 
further pursue the comparison with experimental data in 
this paper leaving that for a future study. 

IV. ANALYTICAL RESULTS IN r s < 1, T/T F < 1 
LIMIT 

We have also carried out an analytic calculation of the 
temperature-dependent 2D effective mass in the leading 
order dynamically screened interaction. This turns out 
to be an extremely difficult task due to the highly compli- 
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FIG. 6: Calculated m* /m at fixed values of temperatures. T 
is in the units of Tf at r 3 = 1. 

cated nonanalytic structure of the integrand in Eq. Q . or 
equivalently Eqs. © or (|T0|> . It is only possible to carry 
out our analytical work in the high-density (r s -C 1), 
low temperature (T/Tp <C 1) limit. It is well established 
that at zero temperature one can do an r s expansion for 
the quasiparticle self-energy in the r s <C 1 limit since 
the leading order contribution in r s comes only from the 
ring diagrams, which are exactly what we calculated in 
our theory (Fig. In this sense, RPA is exact in the 
high-density limit, and one can calculate an exact quasi- 
particle effective mass from the diagrams of Fig. ^in the 
r s — ► limit. In the current work we carry out an exact 
expansion in the r s — > and T/Tp — > limit. At finite 
temperatures, our analysis shows that it is valid to ex- 
pand self-energy and effective mass in r s and T/Tp in the 
r s « 1, T/Tp <C f region. Again we prove that the lead- 
ing order contribution in r s and T/Tp also only comes 
from the ring diagrams, and therefore RPA remains exact 
in this limit. Out calculation shows that 

(13) 

where A(r s ), B(r s ) and C(r s ) are functions independent 
of termpature, and B{r s ) sw Bq > and C(r s ) < for 
r s « 1. 

In Eq. the first term A(r s ) is responsible for the 
nonmonotonic behavior of r s dependence of m* at zero 
temperature, i.e. when T = 0, r s <C 1, effective mass first 
deceases with increasing r s , and then increases. This cor- 
responds to our zero temperature effective mass curve as 
a function of r s in r s <C 1 region, which is shown in Fig.0 
We mention here that this nonmonotonic zero tempera- 
ture r s dependence of the effective mass has already been 
found by previous works£*£*2*S*£. 

The second term in Eq. (|13|l accounts for the leading- 
order temperature correction to the effective mass, which 
is of more interest to us. Our analytical calculation 
shows that B(r s ) w Bq > for r, C 1, ensuring that 



the leading-order temperature correction, in agreement 
with our numerical results, enhances the effective mass 
renormalization in a linear manner as T — > 0. Moreover, 
the linear temperature coefficient is independent of r s as 
r s 1 is also in good agreement with our numerical re- 
sults in this region, as shown in Fig. and discussed in 
section HTll 

The subleading temperature correction, shown as the 
third term in Eq. I|13|l . is negative. This correction, com- 
bining with the leading order linear temperature correc- 
tion, produces the peak in the effective mass temperature 
dependence. This again agrees qualitatively very well 
with our numerical findings in the r s <C 1, T/Tp <C 1 
region as shown in Fig. [S] Of course it is very difficult 
to determine whether the subleading temperature cor- 
rection is T 2 or T 2 lnT dependence just by examining 
the numerical results, but the sign of this correction is 
certainly negative. 



V. CONCLUSION 

Our most important new result is the unanticipated 
non-analytic linear-T enhancement of the quasiparticle 
effective mass at low T/Tp and for all densities. This re- 
sult transcends our specific GW approximation scheme 
since it persists for r s <C 1 where our approximation 
is exact. Since all quantum many-body renormalization 
must vanish in the classical high temperature limit, it fol- 
lows rigorously that m*(T) must be nonmonotonic with 
a peak somewhere at an intermediate temperature as 
shown in Fig. 0] We point out, however, that this non- 
monotonicity would not be easy to observe experimen- 
tally since the quasiparticle is unlikely to be well-defined 
at finite values of T/Tp (~ 0.2 — 0.8) where the peak of 
m*(T) lies. On the other hand, it should be possible to 
experimentally verify our predicted non-analytic linear in 
T enhancement of the quasiparticle effective mass at low 
T/Tp. 

Finally, we comment on the approximations used in 
our calculation. First, our theory leaves out quasi-2D 
form factor (and related solid state physics) effects which 
are straightforward to include^ by appropriately modify- 
ing the bare interaction V q in the theory, and would not 
lead to any qualitative changes in the results (but would 
reduce the magnitude of the mass renormalization by a 
factor of 1.2 to 2 depending on the electron density). Sec- 
ond (and more importantly), our use of the leading-order 
GW-RPA approximation, which is exact only in the high 
density (r s <C 1) limit, is open to question. Although we 
believe that at finite temperatures the GW-RPA approx- 
imation becomes more accurate (and our quasiparticle 
energy calculation of Eq. Q approximately incorporates 
some vertex corrections going beyond the leading order 
expansion in the dynamically screened interaction^*^), 
our principal rationale for carrying out the GW-RPA 
many-body calculation is that (a) it is the only system- 
atic many-body perturbative calculation that is feasible 
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for interacting quantum Coulomb systems; and (b) RPA, 
while being exact only in the weakly interacting r s « 1 
limit, is known to produce qualitatively reasonable re- 
sults even in the strongly interacting (r s > 1) regime, as 
demonstrated by the agreement between RPA and exper- 
iments in 3D metals (r s ss 3 — 5) and in 2D semiconduc- 
tor systems (r s w 1 — 10). The fact that our predicted 
nonanalytic low temperature many-body enhancement 
of effective mass systematically persists to the r s <C 1 
regime shows the generic validity of our results. In ad- 
dition, RPA self-energy calculation should become more 
accurate as r s increases (i.e. decreasing density) for a 
fixed non-zero T because RPA is exact at any density for 



T/T F > 1. 

In this context we emphasize that the RPA self-energy 
calculation (i.e. our effective mass calculation based on 
the diagrams of Fig. ^) is an expansion in the dynamically 
screened Coulomb interaction which becomes equivalent 
to an expansion in r s only in the r s — ► limit. The RPA 
self-energy at arbitrary r s may not be an expansion in 
r s at all, but in some other effective parameters. Even 
in the high-density r s — > limit, the effective expansion 
parameter turn out to be r s /j where 7 is a number of 
order 15(5) in 3(2) dimensional systems. 

This work is supported by NSF-ECS, ONR, DARPA, 
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